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In this work I introduce a powerful analytical method to analyze states of non-equilibrium po¬ 
lariton Bose-Einstein condensates (BEG). It is shown that the condensate wave functions carrying 
dark solitons and half-dark solitons can be expressed in terms of the natural parameters of the non¬ 
equilibrium system, i.e. the growth and decay rates and the interaction parameters. The derivation 
is based on a non-conservative Lagrangian formalism, which provides ordinary differential equations 
for the parameters of the dark soliton solutions in their dynamic environment. Explicit expressions 
for the stationary solution are stated, which significantly differ from those observed in atomic BEG. 
Subsequently the method is extended to spinor condensates, which yields ordinary differential equa¬ 
tions for the parameters of half-dark solitons. Finally a stationary case with explicit expressions for 
half-dark solitons is presented as well. 
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Introduction. - Polariton Bose-Einstein condensates are 
prime examples for non-equilibrium systems of coherent 
quantum matter waves that can be easily controlled ex¬ 
perimentally in contrast to atomic BEG [ij-y . However, 
the theoretical description in the mean-field regime is 
more complicated due to the non-conservative charac¬ 
ter of the underlying partial differential equations, and 
so far effective and simple methods to analyze solutions 
have been lacking. BEG have been studied for decades 
particularly due to their intrinsic nonlinearity and rela¬ 
tively simple mathematical description, particularly be¬ 
cause these states of coherent quantum matter support 
excitations such as solitons 5|-|8|| , quantum vortices, vor¬ 
tex rings or giant vortices which essentially require 

such nonlinearities 1^. It turns out that these excita¬ 
tions depend significantly on the unique features of the 
underlying particles forming the BEG. 


Polaritons are quasiparticles formed within semicon¬ 
ductor microcavities by external illumination via a laser 
beam, which can be created locally at will either co¬ 
herently or incoherently 0, i- Polaritons manifest a 
strong coupling between excitons within the semicon¬ 
ductor created by excitation of the electron into the 
conduction band and cavity photons driven by the ex¬ 
ternal light source. Excitons are electrons coupled to 
oppositely charged holes correspondingly in the semi¬ 
conductor and are held together by the Goulomb force 
between them 0]. The so called bright states of ex¬ 
citons can interact with light fields 0 i El and as 
the light field oscillates rapidly within the microcavity 
strong coupling between excitons and the cavity photon 
can occur and so together they form the so called polari- 
ton quasiparticle within the semiconductor microcavity 
[1, El- Although electrons and holes are fermions, the 
composed polaritons obey half-integer spin and therefore 
can condense to the same ground state macroscopically 


and form a quasi BEG. As their mass is 10® times lighter 
than that of rubidium atoms 0 , the critical temperature 
at which they BE condense is much higher and obser¬ 
vation within the Kelvin ranre have been reported in 
GdTe/GdMgTe micro-cavities 0l2| and more recently in 
flexible polymer-filled microcavities even at room tem¬ 
peratures [l8|. This makes these states very attractive 
for applications as the plethora of papers on this mat¬ 
ter shows E3-E1. An interesting feature of polariton 
condensates is that they can show polarization or inco¬ 
herence between the two internal quasi-spin states [ll.E5|. 
which can stem from the polarization of the light source 
SEl- This splitting between the two spin components 
allows the emergence of intrinsically new excitation of the 
condensate wave function such as half-solitons or half¬ 
vortices 0 El, [13. 


One effective theory for the time evolution of the 
polariton-condensate wave function in the mean-field 
re gim e is the complex Gross-Pitaevskii equation (cGPE) 


_ _ In contrast to the Gross-Pitaevskii equation 

for equilibrium BEG H, H, 31|, the cGPE includes com¬ 
plex spatially dependent terms inducing a local gain (for 
positive sign of the term) or loss (negative sign) of the 
mass of the condensate wave function. This change in 
mass for non-equilibrium condensates stems from the rel¬ 
atively short lifetime of polaritions that is according to 
state-of-the-art experiments between 1 — lOOps depend¬ 


ing on the quality of the microcavity confining the light 
held 0, 0, [sl, [33 . The heeting existence of polaritons 
suggests crucial differences for the possible excitations 
of the condensate, in particular the pattern formation 
of dark solitons in contrast to those observed in atomic 
BEG 0E3- E.g. the gain of mass increases the speed 
of sound and consequently the hrst emergence of exci¬ 
tations in superhuid how [5j. While analytical solutions 
describing the dark soliton have been derived within the 
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classic Gross-Pitaevskii theory [IJ, [3^ a rigorous treat 


merit of dark solitons in non-equilibrium polariton con¬ 
densates will be provided here. Introducing the explicit 
theory for non-equilibrium dark solitons utilizing a non¬ 
conservative Lagrangian formalism and the Kantorovitch 
method [s^ allows fundamental statements about the 
amplitude, width, phase and particularly their interplay 
with pumping and decay parameters. Furthermore I ex¬ 
tend this method to spin sensitive polariton condensates 
and provide the analytical treatment for half-dark soli¬ 
tons, i.e. to half-dark soliton solutions showing strong 
pairwise decoherence as predicted in Q. 

The model.- The spinor polariton field xp = 
is governed by cGLEs coupled to a rate equation for the 
excitonic reservoir hr Emilia, 


idtili± = 


( 1 ) 


+ ai(|'0±P + nR) + 

-ba2|V'Tp + |(7’^fl-r)^V’±- (2) 

ai denotes the self-interaction strength and a 2 = —O.lai 
the interaction strength between components, m is the 
effective mass of the polaritons. 7 denotes the scattering 
rate of the reservoir into the condensate and F the decay 
rate of condensed polaritons. The reservoir dynamics are 

dtUR = P - nR{TR + + \ipz^\^)). (3) 

Here P denotes the pumping distribution and F^; the 
reservoir decay rate. For fast reservoir relaxation Tfi » F 
the reservoir dynamics are much faster than that of the 
condensate. Therefore the reservoir Q can be approxi¬ 
mated in a first step as [ 1 ^ 

P ~ -hydV'iP + IV'tP))- (4) 


Single non-equilibrium condensate. - For coherent spin 
states the set of equations simplifies to 28, 2 ^, 4^ 


^ + a'ld^P + nR)+ 

+ ^{'inR-r)'^ilj. (5) 

Here we approximate the reservoir density by 


of ([5]) simplifying the reservoir and decay of polaritons 


as m 




i.e. 


^ (yrifi - F) ~ /3i|'0p). (7) 

I note that in the case of bright solitons has been 
studied for the cGPE with cubic-quintic nonlinearities 
and I employ and extend this methodology in the follow¬ 
ing. To proceed we set the sign of the linear complex 
term S to be positive and so e = — |e| models the density 
dependent decay. Further we neglect the blueshift from 
the reservoir. Rescaling (l5|) r = y/2m/h ■ x, z = ht and 

Ip —s> —^E and using IT]) yields 
V “1 


id^E = -AE+\E\'^E + i6E + ie\E\^E. ( 8 ) 


with the ID Laplacian AE = d^E and e = /3i/ai. For 
repulsive self-interactions given that 5 = 0 and e = 0 
has a general analytic dark soliton solution [l^, 0, |^, 
43 , i.e. 


E{z, r) = {A - tanh(H(r — Cz)) iD)-eyi];>{iz-\-Kr), (9) 

where H, P, C, D, K are constants. 

Following 3^, 4^ 441 the Lagrangian of our system L 
can be separated into the conservative part Lc (real¬ 
valued energy terms) and a dissipative part inducing 
the non-equilibrium properties of the polariton conden¬ 
sate due to pumping and decay of the polaritons, Lq 
( complex-valued terms). The total Lagrangian is given 
by L = Lc -|- Lq with the components 


Lc = ( ^ (Ed.E* - E*d,E) + |VP|2 + ^ 


( 10 ) 


and 


Lq — i -|- £ 




( 11 ) 


By Hamilton’s principle vanishing variation P(J /(Lc + 
hQ)dzdr) = 0 is obtained for the minimizer, i.e. a solu¬ 
tion of the Euler-Lagran ge e quation (jS]) that corresponds 
to the Lagrangian L via |45l | 


E 


d f d 


L- 


d 


dr] \ dE* dE* 


L =0. 


( 12 ) 


We optimize over the transverse coordinate (Kan¬ 
torovitch method) and get [^, [44 1 


For repulsive condensates the coefficient of the real¬ 
valued nonlinear term is > 0. One can further apply 
the Keeling/Berloff approximation for the complex terms 



2Re 



dr 



dr]dE'* ) 


(13) 
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with r] = r, z. By denoting Lc = drLc the Kan¬ 
torovitch formula (1131) can be rewritten as 

A (AA = 

dz V di^' ) dv 

/ OO 

dr{i5E + ie\E\^E)—. (14) 

-OO 

Here v = A{z), W{z), P{z),C{z) as we are using a dark 
soliton ansatz similar to dH]) for the limiting case U —>■ 0 . 
To show the non-equilibrium properties let us extend the 
analytical expression as follows. 

E{z,r) = 

= A{z) ■ tanh exp{iP{z) + iC{z)r). (15) 

In order that and the r.h.s. of m is integrable 
we renormalize the integral by subtracting the constant 
stationary solution, the procedure for renormalizing the 
nonlinearity of equilibrium solutions is outlined in Hi- 
To obtain the non-equilibrium solutions I suggest by 
analogy to consider the relative pumping strength/decay 
of the dark soliton solution compared to the constant 
ground state (GS). Thus 


= ( 2 AdzE* - E*d,E) + 


\VE\ 


2, 


(16) 


U^^ = i{S{\Ey-A^)+s 


{\Ey-A^y 


(17) 


Using m we obtain the explicit expression for U™ (see 
supplemental material eq. 1). We plug this expres¬ 
sion into the l.h.s of eq. Ildl and obtain four equations four 
the four unknowns v. Variation over A yields employing 
the abbreviation ' = 9, 


To satisfy (EU we set (7 = 0 and so we effectively consider 
a set of three equations for the three unknows A, W, P. 
For vanishing time dependence of W and A, {—A^W)' = 
0. Thus one obtains a simple formula for the amplitude 
of the dark soliton in the non-equilibrium condensate in 
terms of linear pumping S and density dependent decay 
e = — jej parameters, 


A = ± 


V^' 


( 22 ) 


This confirms the intuitive behavior of the polariton con¬ 
densate: As the lifetime of polaritons (for different re¬ 
alizations and e.g. microcavities) decreases the ampli¬ 
tude of the stationary state decreases and as the pump¬ 
ing strength increases the amplitude of the excited state 
does. In addition for the width of the soliton we get using 
dUl), (HH) and (ESI), 


This formula shows that the higher the pumping strength 
i5, the smaller the width of the dark soliton and stronger 
decay implies a wider soliton. Using the formulas for the 
width and the amplitude we derive the time dependent 
phase P given through (IT51) . i.e. 


P{z) 


+ A^W^)z 3 <5 

-^-h const. = —z + const. 

3W^ 4 £ 


(24) 


Non-equilibrium spinor condensate.- Let us now con¬ 
sider the full spinor system (E)) for one spatial dimension 
neglecting the secondary effects of TMTE splitting and 
energy relaxation [^. As in the previous derivation we 
use the Keeling/Berloff approximation (EJ and generalize 
this expression to the spinor condensate case by 

^ {^nn -T)^ + ea'y\i,+ y + |V^_n). (25) 

Applying a similar rescaling as for (El), we obtain 


P' = 


2 -h 2A^W‘^ - 3(72IU2 

3iU2 


and by variation regarding W we obtain 
lU = + 

VA^ - 3(72 _ 3p/ 
Variation regarding P yields 


td,E± = -AE± + {\E±y+^\E^y)E± 

+ i6E±+iei\E±y+j 2 \E^y)E±. (26) 

Here 7 = 02 / 0:1 is the magnitude of the elastic cross¬ 
scattering between the two spin components and 72 mod- 
(19) els misalignment between self- and cross-spin decay. The 
associated non-renormalized Lagrangians (for dissipative 
and conservative part) for the spinor condensate are 


{-A^wy = -2A^SW -h ^A^eW, 

0 

(20) 

f b - Eld^E^) + 1 VUfc | 2 + 

fce{+,-} ^ 

and variation over (7 gives 


+ ''-^\+i\E+y\E.\y (27) 

- AA^CW = 0. 

(21) 

^ / 
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including cross-interaction terms in contrast to the co¬ 
herent states formalism from above, and 

fcef-h,-} ^ ^ 

+ te^2\E+\^\E.\\ (28) 

which via the last term in (E51) takes into account the 
effects of cross-decay or interspecific competition between 
the two spin components. Let us now adapt ansatz m 
to the spinor case and write 'L' = {E+, E-) with 


that for the spinor case we again write Lc = drhc 
and Lq = drhq. Hence, the integral of the conser¬ 
vative Lagrangian is formally 


j^fren _ 


dr ^ 
feG{+,-} 


2 

+ |vi?fcp ■ 


{Ekd.El - E*^d,Ek) 

iw-Air 


iC (30) 


with the approximation, which includes the renormaliza¬ 
tion introduced in 1 ^, i.e. 


E±{z,t) = 

( r i rn \ 

= A±{z) ■ tanh exp{iP±{z) + iC±{z)r). 

(29) 

In this ansatz component -|- is assumed to carry a dark 
soliton with density minimum at r = — rg while compo¬ 
nent — carries a dark soliton with a minimum at r = rg. 
The form of the ansatz resembles the half-dark soliton 
splitting observed e.g. in [1, I further assume both 
solitons to be far away from each other, i.e. rg ^ 0 , to ob¬ 
tain converging integrals particularly for the energy den¬ 
sity of the cross-interaction term 'y\E^\'^\E±\'^. Thus we 
can approximate the cross-component effects via a dark 
soliton solution in one component and a spatially homo¬ 
geneous solution A± in the other component, which is 
justified by the relation tanh^(r) ~ 1 when r ^ 0. Note 

I 


pCO 

C- dr{\E_f - Al)\A+\^ + 

Jo 


+ 


dr{\E+f^ 


(31) 


In analogy the non-conservative part is 


L 


/ren 

Q 





+ iej2C, 


(32) 


which takes an additional cross-decay term e"f 2 C into ac¬ 
count. The result for the integrated Lagrangian for the 
half-dark soliton spinor solution is 


j^fren _ 


2Al{2 + WliAl - 3(C2 -p Pi))) 2Hi(2 -H ITi(Hi - 3(Ci -f Pi))) 


314^4 


31T- 


-jAlAliW+AW-). 


(33) 


For the generalized spinor condensate the parameters of 
variation are iz = A±{z),W±{z), P±{z) and C±{z). By 
variation regarding A± and W± we obtain the corre¬ 
sponding parameter equations 


Pi = 


and 


2 -h 2AlWl - 3ClWl - \i{A\Wl A AlW±W^) 


3Wl 


W± = ± 


2Al - 6Cl - 3^2 7 - 6 P' 


(34) 

(35) 


Variation regarding P yields 


i-AlW±y = -2Al6W± + -AieW± + 

o 

+ j2e2AlA‘^^{W± + W^) (36) 


while variation over C implies 

- 4.AlC±W± = 2e-i2A\A\ log(2)(VF^ - VF|). (37) 

As 7 —>■ 0 and 72 —0 we restore the coherent case dis¬ 
cussed above. For stationary density profiles and assum¬ 
ing the width of both half-dark solitons to be equal (1571) 
implies that C± = 0 as long as A± 7 ^ 0 and W± 7 ^ 0. 
Then by using eq. [301 we obtain the simple solution for 
the coupled amplitudes of half-dark solitons, 


A± = ±1 


35 


e( 2 - £ 72 )' 


(38) 


Eq. [38] shows that the cross-decay reduces the amplitude 
of each component. The corresponding expressions for 
the width and the phase of the non-equilibrium half-dark 
solitons are provided within the supplemental material 
47| . Higher interspecific competition 72 or lowered cross- 
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interactions 7 widen the half-dark solitons within each 
spin component of the polariton BEC. 


Conclusions.- In this work I have introduced a non¬ 
conservative Lagrangian treatment of non-equilibrium 
polariton condensates within semiconductor microcav¬ 
ities. I have discussed the case of dark solitons and 
half-dark solitons and have derived ordinary parameter 
differential equations describing their dynamical behav¬ 
ior. Assuming the system to obey a stationary density 
profile explicit two parameter expressions for the am¬ 
plitude, width and phase of the non-equilibrium dark 
soliton could be stated. These expressions only depen¬ 
dent on the pumping strength S and decay parameter 
e of the non-equilibrium condensate. This significantly 
contrasts the notions of dark solitons in atomic BEC, 
where the amplitude, phase and width relate to the 
self-interaction strength, the chemical potential and the 
mass of the atoms 1^. Furthermore the case of non¬ 
equilibrium spinor condensates has been considered for 
the simple scenario of two spatially separated half-dark 
soliton as excitations of the corresponding spin compo¬ 
nents. Again explicit parameter equations could be de¬ 
rived by means of a novel generalized non-conservative 
Lagrangian method. In the course of this work a renor¬ 
malization procedure for the non-conservative terms has 
been introduced in analogy to the well-established for¬ 
malism for energy renormalization. 
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